Phase-imprint induced domain formations and spin dynamics in spinor condensates 
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We demonstrate that certain domain structures can be created both in ferro- and antiferro- 
magnetic spinor condensates if the initial phase is spatially modulated. Meanwhile, spin dynamics 
of the condensate with modulated phases exhibits exotic features in comparison with those of a con- 
densate with a uniform phase. We expect that these phenomena could be observed experimentally 
using a phase-imprinting method. 



PACS numbers: 03.75.Mn, 03.75.Kk, 75.45.+J 



Spinor Bose-Einstein condensate (BEC) has been at- 
tracting growing attentions in the last decade since it 
displays a variety of exoticphenomena associated with 
its spin degree of freedom [J 0, 0] • Spin-domain forma- 
tion and spin dynamics are definitely among such topics 
of particular interest. 

Early in 1998, soon after the experimental realiza- 
tion of the spinor BEC the MIT group investigated 
the miscibility of different spin domains in the spinor 
23 Na condensate [||. The spin-dependant interaction be- 
tween 23 Na atoms is antiferromagnetic (AFM) and do- 
main structures were pre-created by applying a gradient 
magnetic field. However, the \mp — ±1) domains be- 
come almost-completely miscible as the gradient field is 
turned down, indicating that spin domains is hard to be 
formed spontaneously in the AFM spinor condensate. A 
new experimental result further rules out spontaneous 
domain formation in 23 Na condensate Q. Very recently, 
a theoretical work demonstrates that the homogeneous 
magnetic field can lead to spatial modulational instability 
in AFM condensates, followed by the generation of spin 
domains [5J. In contrast, the spontaneous domain forma- 
tion has been observed in the ferromagnetic (FM) 87 Rb 
condensate using an in-situ phase-contrast imaging Q. 
This is a pioneering approach in exploring spin domains 
in spinor bosons. Theoretically, the domain formation 
is attributed to the FM interaction between 87 Rb atoms 
which leads to spontaneous polarization in the ground 

state 0,1 EES- 

Spin dynamics in spinor BECs, usually referring to 
evolutions of spin populations in different spin compo- 
nents mm mm, 

arises directly from spin exchange 
collisions. Taking the F = 1 manifold for example, the 
collision process can be expressed as |tof = 1) + \mF = 
— 1) <-> 2\mp — 0), which naturally hold the conservation 
of total spins. Owing to the quantum nature of BECs, the 
collision process is coherent so that it leads to oscillations 
of spin populations. Such coherent behaviors have been 
observed experimentally in 87 Rb condensates of both the 
F = 1 [l3[ and 2 [l4| manifolds, and very recently in 
the AFM 23 Na condensate So far, theoretical expla- 



nations for the spin dynamics are mainly based on the 
assumption that each component shares the same spatial 
wave function [l^, [H, [H| > called the single-mode approx- 
imation (SMA). It seems that the SMA works well for the 
AFM BEC H, but becomes invalid for the FM conden- 
sate due to its domain structures. We have ever proposed 
a two-domain model to account for the later case and 
suggest that domain formations inside FM BECs bring 
about significant influence on spin dynamics [15j . 

In this paper, we propose a scheme for generating spin 
domains both in FM and AFM condensates, and dis- 
cuss the exotic spin dynamics caused by domain forma- 
tions. According to this scheme, domain structures can 
be created by the spatially modulated phases, which is 
expected to be realized via phase-imprinting method in 
experiments [l6| . Phase- imprinting, as a versatile tool to 
manipulate BECs, has already been used to create dark 
solitons [13, vortices [U£ [l9[ and vortex rings [2(| in 
scalar or two component condensates. One can further 
expect that it applies to spinor-1 BECs as well. 

We start with the mean-field energy functional for a 
spinor-1 Bose condensate, which is expressed as Q 



E = dr 



r £2 



2m 



(1) 



where ?p a denotes the condensate wave function for the 
atomic BEC in the a-th internal state \wlf — a) and re- 
peated indices are assumed to be summed. The Co and C2 
terms describe contributions of the spin-independent and 
spin-dependent interactions between atoms respectively. 
The spin-dependent interaction could be FM if C2 < 
or AFM if C2 > 0. F is the vector of spin matrices and 
Vext( r ) is t ne external trap potential. Then equations of 
motion for the spinor condensate, derived from Eqn. (JTJ) 
via variational principles, are given by Q 

d 

ih—ip +1 = [H + c 2 (n+i + n - n_i)]^+i + c 2 ?/'o?/'li , 
d 
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ih gl^-^ = [H + C 2 [n-i +n - n +1 )]-0_! + c 2 ipl'<j)* +l ,(2) 



2 



where H = -|^-V 2 + V ext + c (n + i + n Q + n-i) and n a 
represents density of the mp = a atoms. Equations ([2|) 
are just the Gross-Pitaevskii (GP) equations for spinor- 
1 condensates, with the condensate wave function ex- 
pressed as ip a (r,t) — \Jn a {r, i)e l9 =( r '*). 

The GP equations have been intensively employed to 
describe spin dynamics and domain instabilities of spinor 
BECs, wherein the phase 9 acts as a crucial factor since 
it reflects the quantum characteristics of the conden- 
sate 

@, 01, HO, 111 EH- However, 9 could play a more 
important role than recognized previously. In previous 
works, 9 is usually considered to be spatially invariant. 
Hereinafter, we look at the case that the phase 9 can 
be spatially modulated. Such an extension is by no 
means trivial. We will show that intriguing domain struc- 
tures and spin dynamics could be induced by appropriate 
phase- modulations. 

We consider a one-dimensional (ID) BEC as an ex- 
ample, supposing that it is confined in the infinitely 
deep square well potential, approximately corresponding 
to the elongated cigar-shaped BEC in experiments [2l| . 
The initial condensate density profile is set to be 
homogeneous-like except at boundaries where it tends to 
zero (22[, as plotted in Fig. (Ha). We note that the initial 
phase is not homogeneous, but has been modulated as 



9 a (r,0) =a(l-20(r-L/2)) 



where 



0(z) = 



1 x < 
x > 



(3) 



(4) 



is the step function and L is the width of the well. 

The above initial conditions are very similar to the case 
of generating dark solitons in a one-dimensional scalar 
condensate, where the phase modulation is realized via 
the phase-imprinting method [It} ■ The scalar condensate 
contains only one component, corresponding to the case 
of a = 1 in our model. For spin-1 condensate, there are 
three individual components with probably different ini- 
tial phases, as described by Eq. ([3]) with a = ±1,0 (23j . 

For the ID condensate, the coupling constants cq and 
C2 can be estimate from the 3D scattering length [2l| . 
As a simple approximation, c^ D /cq D rs l/V, where I is 
the elongated length of the condensate, equivalent to the 
width of the well in our model, and V is the condensate 
volume. In following calculations, I is set to be the unit 
of the length. And the unit of time is defined as t un it = 
2ml 2 jh. The effective interaction coefficients are also re- 
scaled, e.g., Cq = 2ml 2 Co /ti 2 where cq is the ID coupling 
constant. We note that both t un it and c' are proportional 
to the square of the length. 

First, we investigate a FM condensate, such as 87 Rb. 
According to the experiment data in Ref. @, the 87 Rb 
condensate contains 2.1 x 10 6 atoms and the length of 
the condensate is about 334^m. So f unit « 300s and the 
re-scaled coefficients c' and c' 2 are approximately equal 
to 300 and —1.4 respectively. Then dynamical properties 
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FIG. 1: Spin domain formation in the ferromagnetic spin-1 
condensate. 



of the spinor condensate can be simulated using Eqs. (|2|) 
and ([3]), and the results are shown in Fig. [TJ It can be 
seen that the density profiles of the three components are 
changing in different manners after the evolution begins, 
although they are thoroughly in the same shape initially. 
The formation process of a — ±1 domains are apparently 
demonstrated from Fig. QJb) tofljd). On the right half- 
side of the figures, the +1 component overwhelms the —1 
component, i.e., the +1 domain is formed. Meanwhile the 
— 1 domain appears on the left. The two symmetrically- 
located domains are in accord with the specific choice of 
initial phases as given by Eq. , which implies that the 
domains are closely related to the modulated phases. 

Obviously, both the shape and location of domains are 
varying with the time. A very interesting phenomenon 
occurs during the period from Fig. QJd) to (f), where 
the ±1 domains have been exchanging their positions. 
It is worth noting that the two domains exchange their 
positions periodically with time. This kind of dynamical 
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FIG. 2: Local magnetization density at r = 0.1,0.3 and 0.5. 
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FIG. 3: Evolution of populations of different spin components 
and total spins in a ferromagnetic condensate. 



behavior of the condensate is somewhat similar to the 
soliton-like dynamics in FM spinor condensates studied 
by Zhang et al. recently 24 1 . In Ref. [24[ , the soliton-like 
behaviors are attributed to the exchange interaction C2, 
while here they are resulted from the phase modulation. 
As to one domain, the inside magnetization density M — 
n+ — n_ exhibits perfect oscillation. Figure [2] plots the 
local magnetization density M at three different points 
in the left side of the condensate. M at r = 0.1 and 0.3 
exhibits a good periodicity with almost the same period 
within the time of our simulation. The oscillation period 
is dependent on the dimensionality of the condensate and 
the parameter C2. Nevertheless, at the boundary between 
two domains r = 0.5, the magnetization density is equal 
to and remains unchanging, owing to the symmetry of 
initial density profiles and the symmetric choice of initial 
phases. 

Figure [3] illustrates the evolution of total population 
of each spin component (spin population) and the to- 
tal spins. The spin population smears out detailed in- 
formation of inner structures inside the condensate, but 
it is an important factor concerning spin dynamics of 
spinor condensates and has been intensively studied pre- 
viously [HI, HH, EH, EH • It is already discovered that the 
spin population exhibits a quantum oscillating feature; 
this point is also confirmed in our calculations. Moreover, 
the total spin conserves during the evolution, although 
the local magnetization is allowed. Comparing Figs. [5] 
and [31 one can get an interesting point that the period 
of population oscillations is different from the oscillation 
period of the local magnetization density. This reveals 
that the spin population is not sufficient to characterize 
the whole feature of spin dynamics when the condensate 
has certain inner structures, e.g., domain structures. 

It is important to point out that domain structures 
can appear spontaneously in the FM condensate even if 
the initial phases are not modulated at all. This kind of 
domain formation is attributed to the spontaneous sym- 
metry breaking in the spin space cause by the FM interac- 
tion between bosons 0, H, H, OH and has been confirmed 
experimentally in the 87 Rb condensate ((J. The spon- 
taneous domain formation usually leads to some multi- 
domain structures and domains seem distributed ran- 
domly in the condensate @, [25[. However, the phase- 
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FIG. 4: Domain formation in a ferromagnetic spinor conden- 
sate with phases are modulated according to Eq. ©. (a-e) 
show the distribution of each component at different times, 
(f) shows the evolution of the local magnetization density at 
different positions. 



imprinting induced domain structure is strongly depen- 
dant on the modulation of initial phases. Regular initial 
phases can induce regular domain structures. As indi- 
cated above, one can even produce a simple two-domain 
structure so as to facilitate experimental probing. More- 
over, if the initial phases are not that regular as above, 
an irregular multi-domain structure can be produced. To 
demonstrate this point, we suppose that the initial phases 
are modulated in a more complicated manner, such as 



9 a (r, 0) = sin ((4 + a)nr) 



(5) 



where different components are imprinted with different 
phases. Figs. BJa-e) show evolution of the condensate. 
A multi-domain structure emerges and the condensate 
displays very complicated dynamical behaviors. For ex- 
ample, the local magnetization density shown in Fig. [4](f ) 
is still oscillating, but not periodically. 

Another striking difference between the two cases is 
the time scale of producing domains. The characteris- 
tic time for the spontaneous domain formation strongly 
depends on the spin-dependant interaction (t un it = 
%/ \c2\n), whereas that for phase-imprint induced do- 
mains is mainly determined by the size of the conden- 
sate (i un it = 2ml 2 /H). The spontaneous domain forma- 
tion takes a period of relaxation time of nearly 100ms 
in 87 Rb condensates |6j. In a similar condensate with 
I = 334^m, it takes about 0.12£ unit « 36s to form the two 
domain structure resulted from the phase modulation, as 
can be seen from Fig. [5J In this case, domains appear 
much more slowly. However, the domain formation time 
can be significantly reduced by shortening the conden- 
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spin dynamics measurements show good agreement with 
predictions made on the base of the SMA. Nevertheless, 
spin domains may be induced by applying some driving 
factors. For example, a recent theory predicts that spin 
domains can be generated by applying an external ho- 
mogeneous magnetic field to the AFM condensates [5j. 
So a question arises: whether domain formation could 
be induced if modulating initial phases? To answer this 
question, we need to simulate dynamic behaviors of the 
AFM condensate according to Eqs. and ([3]), with the 
parameter C2 set to be positive. We choose the parameter 
Cg and c' 2 to be 100 and 10 respectively and the obtained 
results are shown in Fig. [5l The evolution process of 
AFM condensates is quite similar to the FM case, except 
that the density profile of each spin component is differ- 
ent from that in the FM condensate. This indicates that 
one can really create spin domains in AFM condensates 
by the phase imprinting method, as does in FM con- 
densates. These simulation results await experimental 
validation, for example, in the AFM 23 Na condensate. 



FIG. 5: Evolution for domain formation in a anti- 
ferromagnetic condensate. Rather similar to the ferromag- 
netic case, the domain structure is mainly dependent on the 
initial phase. 



sate length. When the length is reduced to about 10/im, 
the formation time can be decreased to within 100ms. 
An alternative way to reduce the domain formation time 
is to enrich the phase modulation pattern. As Fig. |4jf) 
shows, domains are well developed at t « 0.01t un it , much 
faster than the two-domain case where t w 0.12t un ; t . In 
some sense, enriching the phase modulation pattern is 
equivalent to reducing characteristic length of the con- 
densate. 

In comparison with the case of FM condensate, domain 
formations in AFM spinor condensates, such as 23 Na, is 
more fascinating. A very recent experiment suggests that 
the spatial domain structure could not be formed spon- 
taneously in the 23 Na condensate Q. Meanwhile, the 



In summary, we have established that the phase- 
imprinting could induce spin domains both in FM and 
AFM spinor condensates. Domains come into being after 
the initial phases are imprinted and the domain structure 
is strongly related to the spatial modulation of phases. 
These characters make the phase-imprint induced do- 
main formation differs from the spontaneous domain for- 
mation in the FM condensate. Even more interestingly, 
phase-imprinting offers the opportunity to study spin do- 
mains and related spin dynamics in the AFM condensate 
where domain structure can not form spontaneously. Our 
results demonstrate that the phase engineering can play 
more important roles in manipulating the quantum fea- 
ture of spinor Bose condensates than previously done. 
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